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Abstract 

In this lecture^ we review several cosmological issues associated with the axion. 
Axion solves the strong CP problem and is a good candidate for the dark matter. 
Limits, which are imposed by the value of isocurvature fluctuations fraction in the 
observed CMBR from WMAP data and the domain wall problem are discussed in 
a supersymmetric flat direction axion model. 



^This lecture was given at Xll'th Baksan International School of Particles and Cosmology, 21-27 
April 2003, Baksan Valley, Kabardino-Balkaria. 



1 Introduction 



At present astrophysical data suggest that the observable part of the Universe contributes 
only a few percents to the total energy density. Most of the content of the Universe is 
in the form of the dark matter (25%) and dark energy (70%). In this lecture we will 
review one of the most viable candidate for the dark matter, which is axion particle,^ 
that originally arose as a solution to the strong CP problem. Axion decay constant must 
be very large compared to the electroweak scale, in order to not confront astrophysical 
data, which yields the so-called "invisible axion" models. Since axion interacts with other 
particles with the strength inversely proportional to its decay constant it turns to be a 
good dark matter candidate. 

It is worth to mention that string theory provides many periodic fields, associated 
with global symmetries spontaneously broken down to some discrete subgroups. Al- 
though, naively the decay constants of string axions are of order of ^^^2 ~ 

lO^^GeV, 

which is too large for a QCD axion decay constant, one might hope to construct a QCD 
axion model out of string theory axions. The rest of these periodic fields may yield a 
solution to the dark energy if they get stuck away from the minimum of their potentials. 

Below we discuss some axion models stressing some important issues which arise if 
the QCD axion does exist and contribute to the dark matter. In section 2 wc briefly 
remind the usual QCD story and the originally proposed "invisible" axion models. In 
the section 3 we discuss supersymmetric axion which corresponds to some flat direction. 
Finally we discuss axion isocurvature fluctuations and the domain wall problem in the 
context of the supersymmetric axion model. The last two issues impose a signiflcant 
constraints on any axion model. We review how one can avoid corresponding difficulties. 



2 The Strong CP Problem and The "Invisible" Ax- 
ion. 



The Lagrangian of QCD, besides ordinary terms, contains the total derivative term or 
^-term: 

where G"''^^ is the dual to the field strength tensor. This term is gauge and Lorentz 
invariant, therefore it should be incorporated in the full QCD Lagrangian. However, one 

^Thcio are two other candidates for the dark matter, which are neutralino or LSP and mirror world 
matter, which we do not discuss here. 
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can check easily that this 6'-term is a total derivative since it can be rewritten as 

= Qd^K^ (2) 

with 

= -^.e^'-'^iAld^Al + ^U^AIA^AI). (3) 

This term being a total derivative does not affect the equation of motion and if one would 
have trivial vacuum structure like that in QED this term would not affect any physics 
at all. However, because of the rather complicated vacuum structure of QCD, i.e. the 
presence of instantons, the integral over this total derivative docs not vanish. Therefore, 
this term contributes to the action of QCD and must be retained. 

First, let us note that this term violates CP symmetry, thus it can give rise to 
neutron electric dipole moment. Current experiments set a strong bound on the value of 
9^, which is 

9 < IQ-^. (4) 



It has to be explained why this number is so unnaturally small. Peccei and Quinn [2] 
proposed a very elegant solution to that problem (although not unique). Namely, consider 
the case when the Lagrangian has an additional U{1)pq symmetry. If this symmetry is 
broken at some energy scale, which we denote by /„, there will be NG-boson generated. 
This NG-boson is called axion^. It is massless and it couples to the gauge bosons and 
to the matter fields with the strength which is inversely proportional to the scale of the 
PQ-symmctry breaking f^. Nonperturbative instanton QCD effects, however, give an 
axion a potential of the form, which is well approximated by 



Vqcd{T) = f!ml{T)[l- cos —]], (5) 




breaking U{1) down to discrete subgroup Zjv, with the mass of the axion depending on 
the quark masses and on the temperature of the quark-gluon plasma. This naturally sets 
the coefficient in front of axion-gluon-gluon term in the QCD Lagrangian to zero, since 
minimum of that potential is at a = 0. Thus, PQ solution solves the strong CP problem. 



At zero temperature the axion mass is given by 

mumdnis ^ 

m„m<i + mdms + rusmu fl fl 



2 rriumdms Aqcd f^ml 

= , , 75 (6) 



where Aqcd is the QCD scale and can be related to the pion mass and decay constant 
by ml — 2(m„ + md)h^^cDl fl- ^ ^QCd ^ fa then the axion turns out to be very hght 

^Neutron electric dipole moment can be computed using current algebra methods [1]. 
'''One of the other two known solutions to the strong CP-problem is what is called Nelson-Barr 
mechanism[3]. Second possibility is to set the mass of u-quark to be zero. 
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comparing to the QCD scale. At high temperature axion is even hghter and its mass is 
temperature dependent [4] 

m„(T) = 0.1m„(T = 0)(^)''. (7) 

There are three historically consequential models of the axion : 

• Weinberg- Wilczek axion 

• KSVZ axion (Kim; Shifman, Vainshtein, Zakharov) [5] 

• ZDFS axion (Zhitnitskii; Dine, Fishier, Srednicki) [6] 

In the first model U{1)pq symmetry is broken at the EW scale. There are two Higgs 
doublets and the axion is ^ 

a = -{v^IuKpo - v^Imxo), (8) 

where 0o ^-nd Xo S'l'c the neutral components of the Higgs doublets. The PQ-symmetry 
is broken at v = ^v^ + v^, where v ~ 250GeV. Such an axion is already ruled out by 
the experiment because it would lead to disastrous loss of energy by various cosmological 
objects via axion emission. The other constraint which is even stronger comes from 
accelerator experiments. 

The second model and the third are so-called "invisible" axion models. In the KSVZ 
model one introduces a complex scalar field $ which couples to the hypothetical quark 
field in the fundamental representation of SU (3)c: 

5/: = ^QrQl + h.c. (9) 

$ is supposed to develop large expectation value /„/ \/2 {fa ^ ^qcd) and axion is defined 
to be 

a = /aArg$. (10) 
The low-energy coupling of the axion to gluons is then given by 

such that ^ — >^ ^ -|- -F- in the lagrangian and one may apply Peccei and Quinn solution. 

In the ZDFS model one starts with the additional scalar SM singlet field S. The 
scale of PQ-symmetry breaking is separated from the EW scale and the axion field is like 
in the Weinberg- Wilczek model but with E field added: 

a = ^{v^Im(f)o - v^Imxo + vj:ImT,o), (12) 
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where V = \Jv'^ + v'^ + v^, ~ if I's is large compared to the Higgs expectation value. 
It can be as large as GUT or even Planck scale. 

Last two models are minimal extensions of the standard model, in which one can 
generate effectively the required coupling of the axion to gluons. At the same time in 
both cases one has to introduce a new scale, which is not defined by the model. This 
scale is essentially time independent and constrained by astrophysical and cosmological 
data to be in the range lO^GeV < /„ < lO^^GeV. In the next section we discuss a model 
where one can go a little further then just introducing a new scale, but rather relating /„ 
to other parameters, such as SUSY breaking scale and the Hubble constant. This would 
have an important consequences, which are discussed in the next two sections. 



3 Supersymmetric Case 

All previous models require some extension of the electroweak model. However, one 
may want to consider some supersymmetric version of the original picture. Such version 
would be, for example, an NMSSM model, which also yields a solution to ji problem in 
the MSSM. One adds additionally at least one singlet chiral field with the most general 
superpotential. For a particular choice of couplings this looks like supersymmetrised 
version of previously mentioned ZDFS axion model. One can generate axion-gluon-gluon 
term and apply PQ solution. However, supersymmetric case is different. Besides light 
axion one has additionally another PQ charged scalar, which is the radial part of a 
complex PQ field. This partner of axion is called saxion. Saxion acquires large mass due 
to supersymmctry breaking and may infiuence evolution of the axion field. There are fiat 
directions associated with different fields and one can imagine the following potential^ 

V{S) = ^-H'\S\'- ml/,\S\' + Vqcd{T), (13) 

where H is the Hubble constant, S — |5'|exp(a/|5'|) is the PQ superfield, \S\ is saxion, 
and a is axion. Such potential naturally appears as a result of lifting of the flat direction 
Shy 

• supersymmctry breaking effects in the early Universe, which leads to the masses of 
order oi [8], 

• nonrenormalizable operators, sTipprcsscd by the Planck scale, the lowest of which 
we take as an example in the potential above, 

• zero curvature term of order of the SUSY breaking scale. 
^This is similar to [7]. 
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This potential has a form of of the " mexican hat" , and corresponding U (1) symmetry 
is broken at the energy scale 

h = \S\^{HM,i)l (14) 

at early time. Note, that this effective scale of U{1) breaking is time-dependent until 
H ~ m3/2, when it gets frozen at /„ ~ (r«3/2A4pi)^^^ ~ lO^-'^GeV. We will discuss simple 
model with one inflaton which decays in the end of inflation quickly reheating the Universe 

up to some temperature Tr. In that case Universe becomes first radiation dominated. 
While it expands, it cools and at temperature T < Aqcd the NG boson, which travels 
along the bottom of "mexican hat" gets the mass niaiT < Aqcd) ~ fir'nT-tr/ fa and starts 
to oscillate. This happens at the value of the Hubble constant given by 

H<}cn={^fH„ (15) 

where Hj is the value of the Hubble constant during inflation. Taking Tr ~ lO^GeV, 
and Hi ~ lO^^GeV one obtains 

Hqcd ~ 10-^GeV. (16) 

The initial amplitude of the axion field is ao ^ fa, which corresponds to j- = 6a ~ 0{1).^ 
Therefore from the kinetic term one can get a rough estimate of the axion energy density 
at this moment 

P ~ ^a/a ~ />' ~ (0.1GeV)^ (17) 
which constitutes a tiny fraction of the total energy density 



/2 

Ptot 



10-'^ (18) 



More accurate estimates, including decays of the axionic strings and domain walls lead to 
a larger value by an order of magnitude (for review, see, for example [9, 10]). Axion is so 
light and weakly interacting that it remains being unthermalized, coherently oscillating 
field until present times contributing to the CDM density. 

Since axion exists in the form of the cold matter its energy density decreases slower 
by a factor of a, then that of radiation, with a being the scale factor in Robertson- 
Walker metric. Because of that, this small fraction grows over time and eventually starts 
to dominate the energy density. In the simple model, which we consider, that restricts 
fa to be somewhere lower than 10^^ — lO^^GeV, which is consistent with the value we 
have obtained above. 



^ There is no reason why initially the phase of the axion should be much different from unity, unless 
one employs anthropic principle. 
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4 Isocurvature Fluctuations and Domain Walls 



First issue which shall be discussed is the limit coming from the observed CMB anisotropy 
[11] on the value of isocurvature fluctuations produced by axion. Isocurvature fluctuations 
appear when there are two or more fields, which generate quantum fluctuations during 
inflation. Fluctuations in the density of one of them may be compensated by some 
other field, so that curvature remains unperturbed. These kind of fluctuations thus are 
called isocurvature fluctuations. In the case of massless (or nearly massless) axion and 
massive saxion isocurvature fluctuations are fluctuations in PQ charge or, in other words, 
fluctuations in the phase 9. The size of primordial spectrum of such fluctuations is set 

by 

- = ^a^, (19) 

where Qa is axion contribution to the total density. This can be understood as follows. 
First, one finds that the equation of motion 

P 

9k + 3H9k + -^9 = (20) 
after rescaling 9k ^ ^ and going to conformal time reads as 

0'^ + {k' - -)9k = 0, (21) 



and has a solution 



which gives at superhorizon scales 



H 
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l^'^' -2^-2^- ^^^^ 



Since power spectrum is defined as 



one can see that the size of the phase fluctuations is if//27r/a. 

The observed size of CMB anisotropy is of order of 10~^, which means that 

Q„-^ < 10-5 (25) 
27r/a 
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at least^. Moreover, the observation tell us that fluctuations are almost pure adiabatic 
which restricts (25) even stronger. Although, it is well known that isocurvature fluctua- 
tions tends to decay into adiabatic ones [13], in the case of the axion the limit above still 
is a considerable underestimate. 

However, even (25) impose stringent constraints. Having Hubble constant of order 
of lO^^GeV requires /„ > lO^^GeV, which is far above the upper limit on imposed 
in order to not overdose the Universe. In the scenario described in the previous section 
this problem is naturally avoided, since fa is time dependent. During inflation its value 
is given by [HjMpiY^'^ ~ lO^^GeV, which has just about right size.^ After inflation the 
value of fa decreases with time as until it freezes at (m3/2Mp/)^/^ ~ lO^^GeV. Such 
scenario seems quite plausible because it solves isocurvature fluctuations problem, and 
one does not overdose the Universe at the same time. 

The other question which shall be discussed is the domain walls problem. The 
problem arises as following. When U{1) is broken down to Z„, the vacuum of the theory 
falls into N degenerate, but disconnected and, thus, distinct regions. As it is well known, 
this leads to topological defects which appears as kink solutions between each of the 
disconnected piece of vacuum manifold. Simple example is the Higgs-like potential of 
the real scalar fleld which possess Z2 discrete symmetry. Namely, this theory has two 
vacua = ±i». If there are two regions in space lying in different vacua, one can easily 
obtain kink solution of the classical equation of motion, with boundary conditions 



It easy to see that there is energy density between boundaries, which has maximum at 
a; = 0, and its characteristic thickness would be inversely proportional to the mass of the 
scalar. This object is a simple example of a domain wall. 

Topologically formed domain walls in case of the exact discrete symmetry are ab- 
solutely stable. Since domain walls are two-dimensional objects their energy density 
behaves as 1/a, while the Universe is expanding. This is much slower than that for mat- 
ter (1/a'^) and radiation (1/a*^), thus, leading to cosmological disaster very fast. Since 
we do not observe our Universe to be domain walls dominated, means that if they were 
ever formed their density either was diluted up to cosmologically safe densities, i.e. less 
then one per horizon, or there are no exact discrete symmetries, so that domain walls 
collapse before they dominate Universe.^ 

Below we will discuss both possibilities to avoid domain walls problem. A simple 
■^For more detailed analysis see [12] 

°If one includes higher order nonrenormalizable terms of the form ' jiin — the PQ symmetry is broken 



at higher values /„ > 10"'^ GeV, which relaxes the constraint from isocurvature fluctuations. 

^One other issue is that domain walls, living long enough, lead to a nongaussianity of CMB, which 
is in direct contradiction with current observation [14] . 




(26) 



pi 
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resolution comes in the case when PQ phase transition occurs before or during inflation. 
Parts of the Universe sitting in different vacua grow into very large regions, such that our 
Universe turns out to be living in one of those. We need to make sure, that reheating of 
the Universe does not restore U{1) at any time after inflation, so that the Universe stays 
in the same vacuum until QCD effects turn on. Provided both conditions are satisfled, 
domain walls which might be formed at T ~ ^qcd are not stable topologically. They do 
not separate different vacua and, thus, will decay. 

For that to work we have to require that 

fa>Hi,Tn, (27) 

which seems plausible in the picture described in the previous section, since at early 
times the value of /„ is 

fa ~ iHjMpi)-2 - lO^^GeV, (28) 

when one takes Hj ~ lO^^GeV and Tr ~ lO^GeV as in the simplest chaotic inflation 
scenario. Therefore, the simple construction described in the Section 2 works well to 
solve the domain wall problems too. 

There is one more case one may want to look at. Similar to (5) we can imagine the 
potential of the form 



MS 
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V{S) ^ '-j^ + H'\S\' - ml/,\S\' + Vqcd{T). (29) 

In that case phase transition occurs late, with 

/a~(m3/2MpO^~10^'GeV. (30) 

This automatically solves isocurvature fluctuations problem, since there are none gener- 
ated during inflation. But one is still confronted with the domain walls problem in that 
case. There is still a possibility to made these walls to collapse. For that one needs to 
break symmetry, adding, for example, to the potential 

5y = /cos(^), (31) 

where ji is some scale. However, if one uses the restriction from neutron electric dipole 
moment from one hand, and require domain walls to collapse fast enough to not became 
cosmologically dangerous from another, one arrives to the following criterium 

lO-^GeV <ii< 10-^GeV (33) 



which numerically gives 
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up to factors of order one. This seems rather narrow window, which is hard to satisfy. 
Ahhough, there might be constructed a model, where these hmits are much broader [15]. 

In the end, we should note that the domain walls problem is safely avoided, if the 
discrete group is simply Z\. This case can be constructed, and domain walls will not be 
topologically stable. Therefore, the issue with domain walls is significantly relaxed. 

Finally, we should mention that there are other cosmological scenarios, i.e. those 
where the Universe is matter dominated soon after infiation [16], due to the moduli os- 
cillations, such that Universe reheats several times. First is due to infiaton decay, second 
and further ones are due to moduli decay. There might also be an issue of parametric res- 
onance effects [17], if one imagine that axion potential is very flat. Such scenario solves 
isocurvature fluctuations problem. However, it poses a domain wall problem, because 
after inflation axion acquires large fluctuations due to the instability, which is similar 
to that in Mathieu equations. This case is safe from the domain walls problem only in 
= 1 models^°. 



5 Conclusion 

In this work we briefly went over axion story and two important issues which are 

• the limits on axion models coming from the bound on the value of isocurvature 
fluctuations and 

• the domain walls problem 

Both issues seem to be very restrictive and provide a good test on the viability of any 
axion model. One can construct reasonable supersymmetric models which may safely 
avoid these problems. The domain walls may be cither diluted if PQ transition occurs 
prior to the inflation, or one can break Zjss symmetry by a tiny bit, making one of the 
vacua more preferable then the others. One can construct models, which are domain 
walls problem free. Isocurvature fluctuations problem is solved in a model with time 
dependent PQ scale. More extended discussion of supersymmetric axion models is given 
in [18]. 



Acknowledgments. 

'N = 1 domain walls are discussed in details in [19] 
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